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Abstract. We report on Keck optical BVRI images and spectroscopy of the companion of the binary millisecond pulsar 
PSR J02 18+4232. A faint bluish (V = 24.2, B - V = 0.25) counterpai't is observed at the pulsar location. Spectra of this 
counterpart reveal Balmer lines which confirm that the companion is a Helium-core white dwarf. We find that the white dwarf 
has a temperature of T^g = 8060 ± 150 K. Unfortunately, the spectra are of insufficient quality to put a strong constraint on the 
surface gravity, although the best fit is for low logg and hence low mass (~ 0.2 Mq), as expected. We compare predicted white 
dwarf cooling ages with the characteristic age of the pulsar and find similar values for white dwarf masses of 0.19 to 0.3 Mq. 
These masses would imply a distance of 2.5 to 4 kpc to the system. The spectroscopic observations also enable us to estimate 
the mass ratio between the white dwarf and the pulsar. We find q = 7.5 + 2.4, which is consistent with the current knowledge 
of white dwarf companions to millisecond pulsars. 

Key words. Pulsars: individual (PSR J0218+4232) - stars: neutron stars - stars: white dwarfs 



1. Introduction 

Millisecond or recycled pulsars are almost always found in 
binaries. Their short (millisecond) periods are believed to be 
caused by the accretion of matter from the binary companion, 
which spins up the neutron star, and, by means poorly under- 
stood, decreases its magnetic field. This accretion will happen 
when the companion, in the course of its evolution, overfills its 
Roche lobe. 

When the companion of the pulsar is a white dwarf, one 
would expect the cooling age of the white dwarf to match 
the age of the pulsar. This because the cooling and spin-down 
clocks start ticking at roughly the same time, when the com- 
panion starts to contract to a white dwarf and the pulsar turns 
on following the cessation of mass transfer. 

The age of a spin-down powered pulsar, with period P and 
period derivative P, is given by 



TPSR - 



(n - l)P 



(1) 



where Pq is the initial period and n is the braking index. The 
latter determines how the rotational energy loss scales with the 
period (P oc P^^"). For magnetic dipole radiation (see Lyne 
& Smith 1998 1 we have n = 3, and if we further assume that 
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Pq P, the pulsar spin-down age is given by T^har - P/C^P), 
the "characteristic age" of the pulsar. The uncertainty in the 
pulsar age is caused by n and Pq, where, for the former, values 
lower than 3 will lead to ages longer than the characteristic 
age, while initial periods close to the current period will lead to 
shorter ages. 

Another possible source of uncertainty in the characteristic 
pulsar age is the transverse velocity of the pulsar, which causes 
an apparent acceleration that contributes to the observed period 
derivative. Ignoring this 'Shklovskii' effect (Shklovskii lT970l 
can lead to significant underestimates of the characteristic age 
(Camilo et al. 1994i. 

The cooling age of the white dwarf can be determined from 
its effective temperature and its mass, using a cooling model. 
With the optical detection of several white dwarf companions 
to millisecond pulsars (e.g., Danziger et al. 1993, Nicastro et 
al. ll995l Lundgren et al. 1996a b), considerable effort has been 
spent in constructing such cooling models (Hansen & Phinney 
ll^^Sal Driebe et al. 1998 Serenelli et al.fMB. It was found 
that the main uncertainty in predicting white dwarf cooling 
ages is the thickness of the Hydrogen envelope surrounding the 
Helium-core. Thick envelopes can support residual Hydrogen 
shell burning, which slows the cooling of the white dwarf, 
while thin layers cannot and thus lead to faster cooling. 

Hydrogen shell burning can become unstable, leading to 
thermonuclear flashes, which may reduce the mass of the enve- 
lope until it cannot sustain Hydrogen fusion any longer It was 
already found by Webbink (il975i that these thermal flashes do 
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not occur for white dwarfs with masses below ~ 0.2 Mq, a 
result that subsequently has been confirmed by Alberts et al. 
lfT^6l. Driebe et al. ( 1998 l, Sarna et al. (2000l, and Althaus et 
al. ( I2OOI I. This limit is due to the fact that flashes are caused 
by unstable Hydrogen burning through the CNO cycle (Driebe 
et al. I1999> which does not occur for stars with masses below 
~ 0.2 Mq. 

Above this limit all authors do find flashes, but the con- 
sequences differ Driebe et al. (1998 1 find that the flashes do 
not reduce the envelope mass sufficiently to suppress residual 
Hydrogen shell burning and hence they predict long cooling 
ages for masses also above 0.2 Mq. In a later paper, however, 
Schonberger et al. ( I2000> suggest this result might be invalid 
for closer binaries, in which the mass of the envelope might 
be reduced due to Roche-lobe overflow during a flash. Sarna 
et al. (1200 It used binary evolution calculations to pursue this 
idea further. On the other hand, the Hydrogen envelope may 
also be reduced by nuclear burning during the flashes. Indeed, 
Serenelli et al. (1200 U and Rohrmann et al. ( I2002> find that for 
masses above ~ 0.18 M©, their models experience one or more 
thermonuclear flashes, leading to thinner Hydrogen envelopes 
that cannot sustain nuclear fusion and hence to white dwarfs 
that cool much faster 

The models are best compared with data. In Table|6l we list 
the three millisecond pulsars with white dwarf companions for 
which the mass and temperature are known, and the cooling 
age can be determined. The first system, PSR J1012-I-5307, has 
a pulsar with a characteristic age of 8.9 Gyr, and a white dwarf 
with a temperature of ~ 8600 K. The other two systems, PSR 
B1855-H09 and PSR J0437-4715, both have ages of 5 Gyr and 
white dwarf companions with T^s ~ 4800 K. 

For the cooUng ages of the white dwarf companions of 
these systems to equal to the characteristic pulsar ages, one can 
readily see that there must a difference in the cooling properties 
of these white dwarfs. The systems PSR B1855-^09 and PSR 
J0437-4715 are very similar to each other, but rather different 
from PSR J1012H-5307, the latter being older, but, paradoxi- 
cally having a hotter white dwarf companion. Given the cooling 
models and the masses of these white dwarfs, it seems likely 
that this discrepancy is due to a difference in thickness of the 
Hydrogen layer, those of PSR B1855-H09 andPSR J0437-4715 
being fliin and fliat of PSR J1012-H5307 being fliick. 

Here we report on optical observations of the white dwarf 
companion of PSR J0218-H4232, which, as we will explain be- 
low, is expected to have a mass that is between the masses 
of PSR J1012+5307 and PSR J0437-4715, and thus, po- 
tentially, may help us gain insight in the cooling of white 
dwarfs. In Sect. 2, we outline our current knowledge of the 
PSR J02 18 -1-4232 system. We describe our Keck observations 
in Sect. 3, from which, in Sect. 4, we determine the white dwarf 
temperature and surface gravity, together with radial-velocity 
measurements. The mass of the white dwarf and the age and 
distance of the system are described in Sect. 5 and we discuss 
and conclude in Sect. 6. 



2. PSR J021 8+4232 

PSR J0218+4232, a 2.3 ms pulsar, was discovered by Navarro 
et al. (1995 1. From radio timing observations, they found that 
it resides in a 2 day circular (e < 2 x 10"^) orbit around a 
low mass object, with, assuming a pulsar mass of 1.4 Mq, a 
minimum mass of 0.17 M©. Using the statistical argument that 
the probability of finding the orbit with an inclination / less 
than i() is given by P(i < /q) = 1 - cos /q, we can conclude that, 
with 90% confidence, the mass of the companion is 0.17 < 
Mwd/Mq < 0.44, the range in which one would expect the 
object to be a white dwarf with a Helium-core. 

From considerations of the evolution of binaries contain- 
ing millisecond pulsars with Helium-core white dwarf com- 
panions, one expects a relation between the mass of a Helium- 
core white dwarf companion to a millisecond pulsar and the 
orbital period (see Joss et al. I1987> . This relation has been 
modelled by Rappaport et al. ( 1995'l, Tauris & Savonije (I1999> 
and others. For the orbital periods of PSR B1855-I-09 and 
PSR J0437-4715 (Table|6j, the relation by Tauris & Savonije 
(I1999> predicts companion masses to within the observed er- 
rors. Though this relation is not expected to be accurate for or- 
bital periods less than 2 days, the predicted white dwarf mass 
of PSR J1012+5307, having an orbital period of 0.605 days, 
is within 1.5cr of the observed result. Since the orbital period 
of PSR J0218H-4232 is in between that of PSR J1012-I-5307 
and that of PSR J0437-4715 (5.471 days) it is expected that 
PSR J0218-I-4232 has a mass in between their masses. In fact, 
Tauris and Savonije ( 1999J predict a mass of 0.216 M©. This is 
near the mass of ~ 0.18 Mq, where Althaus et al. (20011 predict 
the transition between thick and thin Hydrogen envelopes. 

From radio observations, a dispersion measure of 61.25 pc 
cm"^ was found. In the direction of the pulsar, the Taylor & 
Cordes (1993 1 model for the distribution of interstellar free 
electrons predicts a minimum distance of 5.7 kpc for this dis- 
persion measure. The quoted uncertainty is about 25%, but may 
be considerable larger in this direction as there are few pulsars 
with which the model could be caUbrated. Indeed, a new model, 
presented by Cordes & Lazio ( I2002> . predicts a distance of 2.7 
kpc. 

The characteristic age of the pulsar is 0.46 Gyr (NavaiTO et 
al. [T995l . Recent radio observations of PSR J0218H-4232 with 
WSRT/PuMa (B. Stappers, priv. comm.) have revelled a proper 
motion, of = 3.97 + 1.73 mas/yr. Even for the largest 
distances we will need to consider, this has negligible effect on 
the characteristic age, and hence we ignore it below. 



3. Observations 

The location of the PSR J0218+4232 system was observed 
with the 10-m Keck I telescope at Hawaii, using the low- 
resolution imaging spectrometer (LRIS, Oke et al. I1995> . 
During the night of 1995 November 21/22, R and I-band im- 
ages were obtained. The following night single deep B, V and 
R images were taken. In both nights the seeing was mediocre, 
I'.'O-l'.'l, and the conditions were not photometric. Indeed, the 
images from the first night were of so much inferior quality 
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Table 1. Log of the Keck imaging observations. 



Object" 


Time (UT) 


fint (S) 




sec z 


1995 November 23 


B,V,R 






PSRJ0218 


7:24-7:56 


1200,900,600 


1.09- 


-1.12 


1997 January 9 


B,V,R,I 






PSR J0218 


5:47-5:56 


120,60,60,60 




1.09 


PG 0231 


5:18-5:33 


2x(20,10,10,10) 


1.03- 


-1.04 


SA 95 


5:36-5:43 


10,10,10,10 




1.13 


SA 95 


7:35-7:42 


10,10,10,10 




1.08 


PG 104 


15:40-15:53 


2x(20,10,10,10) 


1.19- 


-1.22 


PG 132 


15:56-16:02 


20,10,10,10 




1.15 



" Full names are PSR J02 1 8-t-4232, PG 023 1 +05 1 , 
PG 1047+003 and PG 1323-086. 



Table 2. Log of the Keck spectroscopic observations. 



Object 


Time (UT) 


fint (S) 




secz 


1998 December 17 










WD 0030+444 


5:47 


180 




1.11 


PSR J0218+4232 


6:05-9:16 


4x2700 


1.08- 


-1.26 


WD 0518+333 


9:24 


180 




1.04 


G 191B2B 


9:36 


30 




1.19 


1998 December 18 










WD 0030+444 


5:00 


180 




1.10 


PSR J0218+4232 


5:11-8:24 


4x2700 


1.08- 


-1.17 


WD 0518+333 


8:34 


180 




1.08 


G 191B2B 


8:48 


30 




1.21 



than those taken the other night that we did not use them any 
further. 

On the night of 1997 January 8/9 the conditions were pho- 
tometric, and short exposures of the system were taken in B, 
V, R and I-band for photometric cahbration. The seeing was 
again mediocre and varied between 1'.'0-1'.'2. Four standard 
fields (Landoltll9921 Stetson l200 0i were also observed in the 
same bands. Table[nshows the log of the imaging observations. 

Spectra of the companion of PSR J0218+4232 were ob- 
tained with the same instrument, now at Keck 11, on the nights 
of 1998 December 16/17 and 17/18. The 600 lines mm"' grat- 
ing was used, covering the wavelength range 3500-6000 A at 
1.25 A pix"^ The width of the slit was 1", resulting in a res- 
olution of ~6 A. The seeing varied between O'.'8-l'.'O on both 
nights. 

On both nights that the pulsar companion was observed, the 
white dwarfs WD 0030+444 and WD 0518+333 (Reid 1996) 
were observed by way of radial-velocity standards, and the 
spectrophotometric standard G 191B2B (Bohlin et al. 1995) 
was observed for flux calibration. Exposures of a HgKrXe lamp 
were taken for wavelength calibration. All spectra were taken 
at a position angle close to the parallactic one, to minimize loss 
of light and velocity shifts due to differential refraction. For 
PSR J02 18+4232, the angle was chosen such that there would 
be another nearby star in the slit. The spectroscopic observa- 
tions are tabulated in Table|2l 



3.1. Astrometry 

For astrometric calibration, we selected all 85 stars from the 
USNO-A2.0 Catalog (Monet et al. 1998 1 that overlapped with 
a 10 second R-band image obtained on November 23, 1995. 
Of these stars, 57 were not over exposed and appeared stel- 
lar and unblended. Their centroids were measured and cor- 
rected for instrumental distortion using a bicubic function de- 
termined by J. Cohen (1997, priv. comm.). These points were 
fitted for zero-point position, scale and position angle against 
the USNO-A2.0 positions. After rejecting 4 outliers, that had 
residuals larger than 0'.'6, the rms residuals were 0'.'18 and 0'.'22 
in right ascension and declination respectively, consistent with 
the expected errors for USNO-A2.0 positions. 

This solution was transferred to the 10 minute R-band (and 
the 20 minute B-band) image, using 184 (209) stars that ap- 
peared on both images and were stellar, unsaturated and not 
blended. Again the zero-point, scale and position angle were 
fitted and the final residuals were <0'.'03 (<0'.'06) in both direc- 
tions. Thus, our positions should be on the USNO-A2.0 system 
to within 0'.'03. The USNO-A2.0 system is on the International 
Celestial Reference System (ICRS) to within about U.'l. 

This latter uncertainty is much larger than the measurement 
error for our optical position, and also greatly exceeds the pre- 
cision with which the very precise timing position (NavaiTO et 
al. 119951 B. Stappers, 2003, priv. comm.) can be tied to the 
ICRS. Within this uncertainty, our images show a single object, 
'X' hereafter, at a position coincident with the timing position 
(see Fig.Q. This object is also present on the 1997 images. 

3.2. Photometry 

We used the DAOPHOT 11 package (Stetson fTWTt for the pho- 
tometry. First, all images were bias-subtracted and flat-fielded 
using twilight flats. Next, we determined a point spread func- 
tion (PSF) for each object image from a selection of relatively 
isolated and single stars. We determined instrumental magni- 
tudes by fitting this PSF to all stars. These fitting results were 
then used to remove all stars, except those used in determining 
the PSF. Using aperture photometry we determined the magni- 
tudes of these PSF stars and calculated offsets between the PSF 
magnitude and the aperture magnitude. 

Instrumental magnitudes of the standards stars were deter- 
mined using aperture magnitudes. This was necessary since 
the images of the standard fields were deliberately taken out 
of focus to ensure that the stars would not saturate the detec- 
tor. Using values of Stetson J2000> we determined magnitude 
offsets, colour and extinction coefficients to calibrate the 1997 
images. The rms residuals after the calibration are smaller than 
0.02 magnitudes. 

Finally, magnitude offsets between the 1995 and 1997 im- 
ages were determined to transfer the photometric calibration 
to the 1995 images. We find that for B, V and R respectively 
90%, 76% and 85% of the flux actually reached the detector, 
consistent with the presence of thin ciiTus during the 1995 ob- 
servations. The magnitudes of the pulsar companion and other 
stars in the vicinity are tabulated in Table |3 We estimate the 
final uncertainty in the zero points at about 0.02 magnitudes. 
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Fig. 1. Optical images of the location of PSR J0218+4232 and its companion. The left-hand panel shows the 1200 s B image, the 
middle panel the 600 s R-band image, and the right-hand panel a zoom-in on the R-band image. In the latter, the circle indicates 
the timing position of Navarro et al. ( .1995) . For clarity, it has been drawn with a radius of 1", twice the 95% confidence error 
radius. We find one object within the circle and name it object X. Other stars in this field are also named and magnitudes of these 
stars are given in Table|3] Object Y, clearly visible in R but very faint in B, is only 2"west of the pulsar location. 



Table 3. Astrometry and photometry of the companion of 
PSR J0218H-4232 and stars in the field. The nomenclature of 
the stars is according to Fig.^ The uncertainties for star X, the 
pulsar companion, are as listed. For the other stars, the posi- 
tions are good to within O'.'Ol except for those marked with a 
colon, which are good to within 0'.'02. For the photometry of 
the other stars, the presence of no, a single or a double colon 
indicates uncertainties of cr < 0.025, 0.025 < cr < 0.075 and 
cr > 0.075 mag, respectively. These errors are instrumental, 
i.e., they do not include the zero-point uncertainty of the as- 
trometric tie (0'.'2 in each coordinate) or of 0.02 mag in the 
photometry. The B, V and R magnitudes are determined from 
the 1995 images while the I-band magnitudes are found from 
the 60 s image taken in 1997, hence the larger I-band error for 
object X. 





a (J2000) d 










ID 


2'' 18™ 


+42° 


B 


V 


R 


/ 


X 


061364 


32' 17'.' 40 


24.40 


24.15 


23.86 


23.76 




+a'03 


±^.'03 


±0.04 


+0.04 


±0.06 


±0.25 


A 


07!029 


32'08'.'81 


21.79 


20.35 


19.43 


18.65 


B 


0T;066 


32'0a'79 


23.43 


21.88 


20.86 


19.77 


C 


07353 


31'38'.'67 


20.73 


20.00 


19.55 


19.18 


D 


08':700 


31'25'.'35 


23.09 


21.60 


20.62 


19.67 


E 


05!015 


31'58'.'27 


23.46 


22.55 


21.98 


21.44: 


F 


06^:238 


32'47'.'09 


20.64 


19.02 


17.88 


16.58 


G 


05!868 


32'52'.'90: 


17.23 


16.81 


16.42: 


16.15 


H 


071058 


32'58'.'75 


23.30 


22.69 


22.18 


21.41: 


Y 


06U94 


32' 17'.' 00: 


26.37:: 


24.78:: 


23.67: 


22.19: 



3.3. Spectroscopy 

All frames were bias-corrected and sky-subtracted using stan- 
dard procedures. The spectra were extracted using an optimal 
extraction method similar to that presented by Horne J1986> . 
Each extracted spectrum was wavelength calibrated using the 
arc lamp exposures, with a correction as detailed below. Then, 



each individual spectrum was flux calibrated using the obser- 
vation of G 191B2B of that night. 

The spectra of PSR J0218+4232 were extracted using the 
spatial profile determined from the brighter star that we had 
made sure was on the slit. After extraction, flux and wavelength 
calibration, all four spectra of each night were added and av- 
eraged. The average of the two averaged spectra is shown in 
Fig. 12 The Balmer lines H/3, Hy and H6 are clearly visible, in- 
dicating a Hydrogen atmosphere, in line with the object being 
a Helium-core white dwarf. 

As we hoped to determine a radial- velocity orbit, we took 
particular pains in trying to verify our wavelength calibration 
and velocity scale. In the process, we encountered a number of 
problems, which we detail below. As will become clear in the 
remainder of this Section, we have not solved the problems, and 
hence there is an additional uncertainty in the velocities we de- 
rive. While in all likelihood this is smaller than the rather large 
measurement errors for PSR J0218+4232, we still will treat 
our velocities with caution. We stress, though, that the prob- 
lems we encountered should have no effect on the main results 
of this paper, which rely on the temperature determination from 
the spectra. 

As a first correction to the wavelength calibration, we mea- 
sured the wavelength difference between the observed and 
known wavelength of the Oi /15577 night-sky emission line. 
We found that these differences had a rms of 0.7 A. We cor- 
rected all spectra appropriately. 

Next, we determined radial velocities of the two radial- 
velocity standards for each night by fitting a 9 parameter pro- 
file to each line, Hfi up to H9, separately. This profile consisted 
of a parabolic term, added to the sum of a Lorentzian and a 
Gaussian, both having the same centre, but different widths and 
amplitudes. This profile fitted the lines extremely well, but un- 
fortunately we found that the resulting velocities for the differ- 
ent lines were inconsistent with each other. 

After correcting for the bary centric motion of the Earth, the 
velocity determined from the H/3 Une in both spectra of WD 
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Fig. 2. The spectrum of the companion of PSR J0218+4232. Shown in the left-hand panel is the average of the observed spectra 
when shifted to zero-velocity (bottom curve) and the convolved modelled spectrum with T^ff = 8000 K and log^ = 7.0 cm s"^ 
(top curve), shifted 1 yuJy upwards. Superposed to the latter are the observed B and V broadband fluxes (also offset by 1 fiiy). The 
right-hand panel shows the H/3 to He lines from the observed spectrum with the convolved model superposed. 



0518-1-333 was consistent with the results from Reid J1996i . 
while the other lines gave higher velocities. In general, the Hy, 
H(5 and He lines gave velocities that were approximately 56, 78 
and 80 km s ' too red. The velocities determined from H8 and 
H9 had too large an error to be meaningful. 

For the spectrum of WD 0030+444 taken on the first night, 
we found a similar result, with the H/3 velocity in agreement 
with the known value and the velocities of other lines following 
the pattern described above. The results for the second night, 
however, differed from those of the first night. First, the Hy6 
velocity was slow by 25 km s"' , and second, the other velocities 
followed a different pattern. 

To verify our methods, we used the flux calibration standard 
G 191B2B, which is a hot sub-dwarf and which also has the H/S 
to He lines in its spectra, though they are much less wide than 
those of WD 0518+333 and WD 0030+444. Its radial velocity 
has not been determined as accurately, but we can still use it 
for relative calibration. We determined velocities with the same 
method as described above, and found that these followed the 
same pattern as the observations of WD 0518+333 and the first 
observation of WD 0030+444. 

On closer inspection, we found that the observation of WD 
0030+444 on the second night differed in flux from the ob- 
servation on the first night, the second spectrum having 40% 
fewer counts. Using the filter curve from Bessell (J_990j, we 
calculated synthetic B-band magnitudes for both spectra. We 
found that the first spectrum gave a magnitude consistent with 
the literature value, while the second spectrum obviously did 
not. Since there was no difference in airmass and conditions, 
we believe that WD 0030+444 was incorrectly placed on the 
slit, thus causing the inconsistent velocities and the flux dif- 
ference. Therefore, we decided not use the observation of WD 
0030+444 on the second night any further. 

We tried to find a reason for the trend in velocity found in 
the different Balmer lines, but could not find any. Therefore, 
we simply used the velocities found from H/?, Hy, HS and He 
(H8 and H9 had too large an error) to calculated a new em- 



Table 4. Radial velocity measurements of the radial velocity 
standards WD 0030+444 and WD 0518+333 and the spec- 
trophotometric standard G 191B2B. 



Date (UT) 




V (km s ' ) 






WD 0030+444 


WD 0518+333 


G191B2B 


1998 Dec. 17 


78.8 ±3.2 


34.7 + 2.3 


22.3 ± 4.0 


1998 Dec. 18 


35.1+4.4" 


34.7 + 2.6 


9.9 + 4.3 



We believe this radial velocity measurement is in error due to mis- 
placement of this object onto the slit. 



pirical wavelength scale, depending linearly on the wavelength 
calibration found from the arc lamp measurements. This cor- 
rection was applied to all spectra, after applying the shift for 
the wavelength difference found from the night-sky line. 

With both corrections in place, the resulting velocities (see 
Table |3 have root-mean-square errors of approximately 14 
km s"'. We take this as an estimate of the additional uncer- 
tainty in any velocity measurement, but stress that since we do 
not know the cause of the problem with the wavelength caU- 
bration, we cannot be certain it applies to all observations. As 
noted, however, it should not influence our main results, which 
rely on the temperature determination. 

4. Temperature, gravity and radial velocities 

We use the photometry and the spectra of the companion of 
PSR J0218+4232 to constrain its effective temperature and 
surface gravity by comparing them with model atmospheres. 
Furthermore, we determine radial velocities of the two spectra 
and use these to estimate the mass ratio between the pulsar and 
the white dwarf and the systemic velocity of the system. 

We take into account the reddening in the direction of the 
pulsar Schlegel et al HMW) give Eg-v = 0.07 + 0.01 and we 
correct our magnitudes accordingly. This reddening value was 
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also used to correct the observed spectra using relations given 
by O'Donnell ( fT994l . 

The two separate average spectra of the companion, taken 
on the two subsequent days, were compared with model spec- 
tra for pure Hydrogen atmospheres, kindly provided by P. 
Bergeron (2002, priv. comm.), to infer the effective temperature 
and the surface gravity. The grid of models used spanned Teff = 
7000 (500) 9500 K and logg = 6.5 (0.5) 8.0 cm s ^, where the 
number in parentheses is the step between the models. We con- 
volved each model with a cut-off Gaussian to mimic the effects 
of seeing and the slit. This convolved model was then fitted to 
each observed spectrum, using data-points around H/3 up to He, 
leaving the radial velocity v, flux normalisation A and slope a 
to be fitted by minimising the x' merit function, defined by 

^2 ^ ^ ^ fi -/[(!+ ^/cUi] X A(^,74340 A)" j' 

Here fi and cr, are the observed flux and eiTor at wavelength 
Ai, and f[A] is the model flux at wavelength A. We introduce 
the slope a to account for possible errors in the (relative) flux 
calibration. 

For each of the two observed spectra a surface of val- 
ues in the Teff , log g plane, spanned by the modelled spectra, 
was created. For each model spectrum the value of the first 
spectrum was added to that of the second to create a third sur- 
face. In Fig. 13 we show contours and the location of the min- 
imum in this latter surface; the surfaces found from the indi- 
vidual spectra are similar in shape. At the minimum, we find 
v^. = 778.8. As we have 606 degrees of freedom, this in- 

/vnun ^ ' 

dicates that the fit is satisfactory, but that likely we underes- 
timated our errors slightly. The errors in temperature and sur- 
face gravity were rescaled accordingly. The best-fit power-law 
slopes for the two spectra are small, -0. 11 +0. 16 and 0.21 ±0. 19 
for the first and second night, respectively, indicating that our 
relative flux calibration was reasonably good. Indeed, also the 
absolute calibration is fair: the B-band magnitude, using the fil- 
ter curve by Bessell <1990t . from our observed spectra is 24.43, 
within the error of the photometric measurement. 

From Fig.|3j one sees that our spectra constrain the temper- 
ature fairly well (Tgs - 8060 + 150 K), but provide only a very 
weak constraint on the gravity (logg = 6.9 + 0.7 cm s"^). This 
is because the gravity mostly affects the higher Balmer lines, at 
wavelengths shortwards of 4000 A, where our spectra are very 
noisy. 

The photometric colours of the companion can also be used 
to infer a temperature. From atmosphere models by Bergeron 
et al. n^95^, we find a temperature of T^ff = 8560 ± 700 K 
for a pure Hydrogen atmosphere with logg - 8. Serenelli 
et al. J2001> and Rohrmann et al. ( I2002t calculated evolu- 
tionary tracks of the formation and cooling of Helium white 
dwarfs. They also give colours based on model atmospheres 
of Rohrmann J2001> . From these we infer a temperature of 
Teff = 8400 + 500 K over a range of surface gravities (log g = 
6.0-7.5). Both these temperatures are consistent with what we 
found from our spectra. 

The radial velocities that we found from fitting the spectra 
to the models by Bergeron are given in Table|5j These velocities 
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Fig. 3. Confidence intervals on effective temperature Teft and 
surface gravity logg as inferred from our spectra. The plus 
signs indicate the parameter combinations of the models in 
our grid. At each point, the best-fit x^ was determined, leav- 
ing velocity, flux amplitude, and spectral slope free to vary. 
The results were interpolated to create a x^ surface, which 
has a minimum at the position shown with a filled circle, at 
Teff = 8060 + 150 K and logg = 6.9 + 0.7 cm s^^. The three 
contours, at;^'^ plus 1, 4 and 9, delineate the Icr, 2cr and 3cr 
eiTor regions. 

Table 5. Radial velocity measurements of the companion of 
PSR J0218-I-4232. The uncertainties listed are measurement er- 
rors, and do not take into account the possible systematic er- 
rors associated with the problems in the wavelength calibration 
(Sect.E3. 



Date & Time (UT) 


MJDbar 


^orb 


V (km s ') 


1998 Dec. 17 07:45 


51164.3281 


0.544 


-211 ±64 


1998 Dec. 18 06:52 


51165.2891 


0.018 


101 ±78 



" Using the ephemeris of Navarro et al. 1 1995 1: 
To^MlDbar 49150.6086, ^=2.0288461 days. 



are corrected for the barycentric motion of the earth and are for 
the best fit of T^s and logg. The error estimates include the 
extra uncertainty in each velocity due to the uncertainties in 
the temperature, gravity, normalisation, and slope. 

The two radial velocities were fitted to a circular orbit using 
the ephemeris of Navarro et al. ( .1995j . We find a systemic ve- 
locity of y = -57+52 km s ' and a radial- velocity amplitude of 
^WD - 159 + 50 km s"' . Naturally, this is a perfect fit, because 
we use two data points to determine two unknowns, i.e., we 
have no degrees of freedom and thus no information about the 
correctness of the fit. This makes the systemic velocity and the 
velocity amplitude very dependent on any error in the determi- 
nations of the radial velocities. We note, though, that possible 
systematic eiTors in the radial-velocity determinations have less 
effect on the velocity amplitude than on the systemic velocity, 
because, to first order, such errors should cancel out when de- 
termining the velocity amplitude, while being compounded in 
the systemic velocity. 
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From the radial velocity amplitude of the white dwarf and 
the pulsar, 21.3243 + 0.0001 km s"', (Navarro et al.[T995}, it 
foUows that q = Mpsr/Mwd =7.5+2.4. 

5. Mass, age and distance 

With our effective temperature and estimate of the surface grav- 
ity, we can use evolutionary models to estimate the mass and 
age of the white dwarf as well as the distance of the system. 

5.1. White dwarf mass 

The weak constraint that our spectroscopic observations place 
on the surface gravity translate into a similarly weak constraint 
on the mass. Using the models of Rohrmann et al. (p002i, 
which span the range of 0.148-0.406 M©, we find a white 
dwarf mass of Mwd = 0.21+° q4 Mq. 

Similar values are found using the white dwarf evolution 
models of Driebe et al. J1998> . We stress, though, the uncer- 
tainty: the upper 2cr error on the surface gravity is log g = 8 
(see Fig. |3}, which corresponds to a mass of about 0.6 
which is above the range expected for a Helium-core white 
dwarf. 

We can also use the mass ratio, determined above, to esti- 
mate the mass of the white dwarf. If we assume a pulsar mass 
of 1.48 + 0.08 Mq, the weighed average of pulsar masses from 
Tabled (excluding that of PSR J0218-H4232), this results in a 
white dwarf mass of Mwd - 0.20 + 0.06 M©, in line with what 
we determine from the surface gravity. The pulsar mass, as pre- 
dicted by our observed values, is MpsR = l-6^0 7 ^o- 

5.2. Distance and age 

To compare the white dwarf cooling age with the characteristic 
age of the pulsar we need to know the white dwarf's mass and 
temperature, as well as the distance to the system. The latter 
can be found from comparison with the white dwarf's absolute 
magnitude and our photometry, but for this we, again, need to 
know the mass of the white dwarf. 

Since our constraint on the white dwarf mass is weak, we 
invert the procedure. For a given white dwarf model and a given 
mass, we compare our temperature with the predictions from 
the model to infer a cooling age and use the luminosity or the 
absolute visual magnitude to find a distance. 

By comparison of our temperature with the cooling models 
of Hansen & Phinney (I1998a> and the evolutionary models by 
Driebe et al. (1998"l and Rohrmann et al. ("20021 we estimate the 
cooling age and the distance of the white dwarf as a function 
of its mass: we show the result in Fig.|4] Each point represents 
a model with a certain mass. The errors in the distance and the 
age are caused by the uncertainty in the effective temperature 
of the white dwarf and, to a lesser extent, the uncertainty in the 
observed magnitudes. The errors are depicted by diagonal bars 
instead of horizontal and vertical ones, because the uncertainty 
in the distance is correlated with the uncertainty in the age, as 
both depend on the white dwarf temperature. The horizontal 
line depicts the characteristic age of the pulsar. 
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Fig. 4. The white dwarf cooling age as a function of distance. 
For each model, we determine the distance and cooling age for 
our observed temperature and magnitude. Since this, through 
the model, creates a correlation between the distance and cool- 
ing age, we use diagonal error bars. The two Hansen and 
Phinney models have larger errors, because we had to use our 
more uncertain I-band magnitude. The dots are the Rohrmann 
models, the squares the Driebe models and the triangles the 
Hansen and Phinney models. The solid curve denotes the char- 
acteristic pulsar age. 

Figure|4]shows that, for this temperature and observed mag- 
nitude, there is a clear distinction between the models predict- 
ing thin Hydrogen envelopes and those predicting thick en- 
velopes. All models by Driebe et al. J1998> and the two lowest 
mass models by Rohrmann et al. ('2002'i have thick Hydrogen 
envelopes and predict white dwarf cooling ages that are much 
older than the characteristic age of the pulsar. The Hansen & 
Phinney ('1998ai models and the Rohrmann et al. ^2002^ mod- 
els with masses above 0.18 M© have thin Hydrogen envelopes 
and hence predict younger white dwarfs at the temperature 
found for the companion of PSR J02 18+4232. If the age of the 
pulsar is indeed roughly equal to the characteristic age, then the 
white dwarf companion must have a thin Hydrogen envelope. 

The second thing to note is the close agreement with which 
the diff'erent groups predict the distance for a given mass. Since 
the luminosity we observe depends on the temperature, dis- 
tance and radius of the white dwarf, similar distances at a given 
temperature indicate similar radii. Apparently, the radius is set 
by the degenerate core, the thickness of the Hydrogen envelope 
having little influence at this temperature. 

Returning to the comparison between the characteristic age 
of the pulsar and the cooling age of the white dwarf we can 
conclude that Helium-core white dwarfs with thin Hydrogen 
envelopes predict ages that are most similar, but a httle older, 
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to the characteristic age of the pulsar. The best agreement be- 
tween both ages is found for white dwarfs with masses between 
0.19 to 0.30 Mq, corresponding to distances of 2.5 to 4 kpc and 
cooling ages around 0.7 Gyr. 

6. Discussion and conclusions 

We have identified the companion of PSR J02 1 8 +4232 and pre- 
sented photometry and spectroscopy. Our spectra show that the 
companion has a Hydrogen atmosphere and from comparison 
with Hydrogen atmosphere models we find that the object has 
temperature of Tf.fi = 8060 + 150 K and a surface gravity of 
logg = 6.9 + 0.7 cm s"'. 

This temperature and gravity would imply a mass of 
■^WD - 0-2 04 Mq, in line with what one would expect for a 
Helium-core white dwarf. Similar white dwarf masses, 0. 19 to 
0.30 Mq, are found when comparing the predicted cooling age 
of the white dwarf with the characteristic age of the pulsar. 

Furthermore, the comparison between the cooling age and 
the characteristic age shows that this white dwarf companion is 
likely to have a thin Hydrogen envelopes, in agreement with the 
prediction of Althaus et al. C2001j . that the transition between 
thick and thin Hydrogen envelopes is situated around a white 
dwarf mass of ~ 0.18 Mq. 

Thin Hydrogen envelopes are also required to match the 
white dwarf cooling ages with characteristic pulsar ages in 
the systems PSR B1855-H09 and PSR J0437-4715 (Table |6j. 
For PSR J1012+5307, several authors (Alberts et al. 1996i 
Driebe et al. ll998l Althaus et al. 2001 1 have argued that a thick 
Hydrogen layer is needed to explain its old age and high tem- 
perature. Our mass determination is to weak to put a direct 
observational constraint on the location of the transition. On 
the basis of the relation between orbital period and mass, how- 
ever, we expect that the mass of PSR J02 1 8 +4232 is in between 
those of PSR J1012+5307 and PSR J0437-4715. We can con- 
clude that for whatever mass corresponds to an orbital period 
of 2 days, thin Hydrogen layers are needed. 

This empirical constraint may be improved from observa- 
tions of the white dwarf companions of millisecond pulsars 
PSR J0034-0534 and PSR J0613-0200, as these systems have 
orbital periods, 1 .6 and 1 . 1 d respectively, in between those of 
PSR J0218+4232 and PSR J1012+5307. 

The companion of PSR J0034-0534 has been detected by 
Lundgren et al. ( i[T996b |l who place an upper limit on the tem- 
perature of 3800 K. This limit depends on the distance, which 
was inferred from the dispersion measure. Nevertheless, with 
a characteristic pulsar age of 5.9 Gyr, this white dwarf would 
likely need a thin Hydrogen envelope to have cooled to such a 
low temperature. 

So far, the only observed system for which a thick 
Hydrogen envelope is required, is PSR J 1012+5307. We should 
keep in mind, however, the possibility that this pulsar is actu- 
ally very young, i.e., that the assumption Pq P does not hold. 
If this were the case then perhaps no white dwarf would need 
thick Hydrogen layers. 

This can be verified observationally by measuring the tem- 
perature of the white dwarf companion of PSR J075 1 + 1807. 
This binary has an orbital period of only 6.3 hours, shorter than 



that of PSR J1012+5307, and hence the white dwarf should 
have even lower mass. Optical observations by Lundgren et 
al. (1996a) have not detected an optical counterpart down to 
V - 23.5, indicating that the white dwarf companion is rather 
cool {T <: 7000 K for M < 0.2 Mq and J ^ 2 kpc). It would 
need a thin Hydrogen layer to have cooled to such a low tem- 
perature within the characteristic age of the pulsar (8 Gyr). 
Identification of this counterpart and a measurement of its tem- 
perature could decide whether it has a thin or a thick Hydrogen 
layer. 
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Table 6. White dwarf companions to millisecond pulsars. The mass and temperature of the white dwarf in these systems is known 
accurately. Listed here are the orbital period, the mass of the white dwarf and of the pulsar, the temperature of the white dwarf 
and the characteristic age of the pulsar. 

Pulsar Name Pb (i) ^?wd (Mq) Mpsr (Mq) Teff (K) Tchar (Gyr) References 
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PSRJ0218+4232 2.029 O-^i-lol i-^ll'j 8060 ± 150 ~ 0.46 9''^J(rJ_ 

References: (i) Lange et al. <2001> . (2) van Kerkwiik et al. <1996> . (3) Callanan et al. 1 1998 1, (4) van Straten et al. (2001 

(5) Danziger et al. {TW^ . (6) Kaspi et al. {TW^ . (7) Nice et al. (8) van Kerkwijk et al. CTUl , (9) Navarro et al. IT9951 
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Timing observations. 
Mass determination. 
Temperature determination. 



O'Donnell, J. E. 1994, ApJ 422, 158 

Oke, J. B., Cohen, J. G., Carr, M., et al. 1995, PASP 107, 375 
Rappaport, S., Podsiadlowski, Ph., Joss, P. C., Di Stefano, R., Han, Z. 

1995, MNRAS 273, 731 
Reid, 1. N. 1996, AJ 111, 2000 
Rohrmann, R. D. 2001, MNRAS 323, 699 

Rohrmann, R. D., Serenelli, A. M., Althaus, L. G., Benvenuto, O. G. 

2002, MNRAS, 335, 499 
Sama, M. J., Ergma, E., Antipova, J. 2000, MNRAS 316, 84 
Sarna, M. J., Ergma, E., Gerskevits, J. 2001, AN 322, 405 
Schlegel, D. J., Finkbeiner, D. R, Davis, M. 1998, ApJ 500, 525 
Schonberger, D., Driebe, T, Blocker, T. 2000, A&A 356, 929 
Serenelli, A. M., Althaus, L. G., Rohrmann, R. D., Benvenuto, O. G. 

2001, MNRAS 325, 607 
Shklovskii, I. S. 1970, Soviet Astron., 13, 562 
Stetson, R B. 1987, PASP 99, 191 
Stetson, P B. 2000, PASP 112, 925 
Tauris, T. M. & Savonije, G.J. 1999, A&A 350, 928 
Taylor, J. H. & Cordes, J. M. 1993 ApJ 411, 674 
Thorsett, S. E. & Chakrabarty, D. 1999, ApJ 512, 288 
van Kerkwijk, M. H., Bergeron, P, Kulkami, S. R. 1996, ApJ 467, 

L89 

van Kerkwijk, M. H., Bell, J. R, Kaspi, V. M., Kulkarni, S. R. 2000 
ApJ 530, L37 

van Straten, W., Bailes, M., Britton, M. et al. 2001, Nature 412, 158 
Webbink, R. R 1975 MNRAS 171, 555 



